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1. Introduction 



We consider pure SU(2) Yang Mills theory at finite temperature in the confinement phase 
T < T c , which is a crude approximation of Quantum Chromodynamics (QCD). We study a model 
of dyons, in particular the quark antiquark free energy, with the aim to better understand the phe- 
nomenon of confinement. In this work we mainly focus on technical difficulties associated with the 
long-range nature of the dyon potentials. Physical aspects and conclusions are discussed in detail 
in [jlj ^] and another talk given at this conference [Bj. 



2. The non-interacting dyon model 



In Yang Mills theory observables O are given by the path integral 

(°) = ^J @ A 0[A]exp(-5YM[A]). (2.1) 

Since there are currently no methods to solve this path integral for low-energy observables analyti- 
cally, one either resorts to numerical lattice gauge theory or to certain simplifying approximations. 

One such approach particularly useful to obtain a qualitative understanding of certain phe- 
nomena of Yang-Mills theory and QCD is the semi-classical approximation. One expands the 
path integral around classical solutions of the Yang-Mills field equations, for which the action is 
locally minimized, and which are expected to dominate the path integral. A specific kind of a 
semi-classical model and its capability to generate confinement is based on dyons. Dyons are lo- 
calized objects carrying electric charge as well as magnetic charge and are named after particles 



with similar properties [M. The path integration (2. 1 ) is transformed from field coordinates to dyon 



collective coordinates and quantum fluctuations, where a Jacobian emerges. Part of this Jacobian is 
the determinant of the so-called moduli space metric. This moduli space metric has been calculated 
analytically for calorons [|h, following [^[ [71 which are pairs of different-kind dyons. A pro- 
posal for a metric of an arbitrary number of same-kind and different kind dyons was made in 
Numerical investigations of this generalized metric however, indicated certain shortcomings, in 
particular its non-positive-definiteness, casting severe doubts on its usefulness. Our main interest 
in this work is to test a numerical method able to treat the dyon long-range potentials in a proper 
way. Therefore, we study a much simpler model of dyons without any interactions, i.e. where the 
moduli space metric is ignored. 

The key observable we are studying is the free energy of a static quark antiquark pair at sepa- 
ration d = |r — r'l, which is given by 

F Q Q(d) = -T log (P(r)i»V)>. (2.2) 

The Polyakov loop correlator (P(r)P^(r')) is obtained as a statistical average in dyon ensembles 
characterized by the spatial dyon density p and the temperature T (we always consider maximally 



non-trivial holonomy, which seems to be intimately connected to the confinement phase [10]). 
The dyon ensembles we study are neutral, i.e. there is an identical number of dyons and antidyons. 
Individual dyon configurations of these ensembles are given by the randomly and uniformly chosen 
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dyon positions {r ; }. One can show that in such ensembles a Polyakov loop is given by 

P(r) = -sin(-^3>(r) ) , (2.3) 



2T 

where <I> is the superposition of the 0-component of the dyon gauge field in the Abelian limit 

<I> ( r ) = I^S- (2-4) 
j \ r i- r \ 

For a more detailed discussion of these equations and the non-interacting dyon model in general 
we refer to [p], W. 

In the following sections we are mainly concerned with evaluating (2.4) numerically, which 
contains an infinite sum over 1/r long-range potentials. 



3. Long-range dyon potentials and finite volume effects 

The 1/r long-range nature of the dyon potential causes severe problems for numerical simula- 
tions. One expects that rather large volumes are needed, to render finite volume effects negligible. 
This in turn amounts to a huge number of dyons, which is proportional to the required computa- 
tional resources. 

A first attempt to simulate dyon ensembles numerically, which suffers from the just mentioned 
problem, is described in [JTJ . There no dyons are considered in a cubic spatial volume of length L. 
Observables are then evaluated in a cubic spatial volume of length I < L located at the center of 
the larger volume (cf. Fig. |l|a). This straightforward method, however, has certain shortcomings: 
(A) reducing finite volume effects to a moderate level requires i <C L, which drastically reduces 
the volume, in which observables can be evaluated; this clearly increases statistical errors; (B) 
an extrapolation to infinite volume is technically difficult, since it has to be done with respect 
to two parameters (£ and L); (C) when attractive and repulsive forces between dyons are taken 
into account, dyons tend to accumulate near the boundary of the large volume, i.e. translational 
invariance is broken severely. 

A better method for treating long range dyon ensembles, which solves or eases the problems 
mentioned above, seems to mimic infinite volume by implementing periodic boundaiy conditions. 
One considers a cubic spatial volume of length L filled with rip dyons (dyon density p = nj)/L 3 ). 
This volume is periodically repeated in all three spatial directions (cf. Fig. |l|b). With this setting 
one expects (A) that finite volume effects are significantly reduced, since the original volume is 
"surrounded by infinitely many dyons" in every spatial direction; (B) an extrapolation to infinite 
volume only has to be done with respect to one parameter, L; (C) interacting dyons will not accu- 
mulate near the boundary of the volume, because periodicity implies exact translational invariance. 



4. Ewald's method 



To implement periodic boundaiy conditions for our long-range dyon potentials we resort to a 
method proposed first in the context of condensed matter physics [11] ("Ewald's method") and is 
widely used in plasma physics, as well. 
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Figure 1: a) Dyons in a cubic volume of length L (blue), evaluating observables in a cubic volume of length 
I < L (red), b) Dyons in a volume of length L with periodic boundary conditions. 



The central quantity, when one is e.g. interested in Polyakov loop averages or the quark anti- 
quark free energy, is the potential of no dyons in a periodic cubic volume of length L, 

<f(r)=II, q ~ i — ] (4-1) 



r — ry + nL 



(cf. also section ^). Due to the long-range nature of the individual dyon potentials, J2nez 3 ■ ■ ■ con ~ 
verges far too slowly for any straightforward efficient numerical evaluation. Ewald's method solves 
this problem by splitting the sum into a short-range part and a long-range part, 

<D(r) = <£ Short (r) + <I> Long (r) (4.2) 
<D Sh °«(r) =Vy -. q -i r erfc ( If^+g^ (4 . 3 ) 

OLong (r) = f£ £ g 1J_ e *(r-Tj) e -AV/2_ (4 4) 

For the long-range part the sum is over momenta k = (27r/L)m with m G Z . The arbitrary param- 
eter A controls the trade-off between the short-range and the long-range part. Both parts converge 
exponentially fast, due to erfc(|r — r ; - + nLl / \/2X) and e~ x ~ k ~l 2 , respectively. A detailed and ped- 



agogical derivation of the splitting into short-range and a long-range part can be found in [ 12]. 



4.1 Evaluating the short-range part 

To determine Polyakov loop averages, one typically computes the dyon potential at M <=< V 
points r distributed throughout the volume. Due to the exponential suppression by 
erfc ( | r — r 7 • + nL | / y/2X ) , it is sufficient to consider dyons inside a sphere with radiu s r max <=c X and 
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center r, i.e. 



* Sh ° rt ( r )= I r J y ,„ r , afc f |r "%t nL| V ^ 



nJ,|r-r ; -+nL|<r lr 



r ; - + nL| I y/2X 



Consequently, the total computational cost is 0(VX 3 ). 

In practice X 3 <C V is chosen. For an algorithm scaling according to &(VX 3 ) one clearly 
needs to determine all dyons close to a point r without iterating over all rip oc y dyons. To this end, 
one divides the original volume into small cubic subvolumes. Technically, the volume is a list of 
subvolumes, where each subvolume corresponds to another individual list of those dyons located 
inside. For any point r one can then easily determine all subvolumes inside or at the boundary 
of the corresponding sphere with radius r max and then iterate over the relevant 6(X 3 ) dyons. The 
approximate sphere of subvolume copies ("LEGO® sphere") is shown in Fig. ^[ Note that the 
radius of the LEGO® sphere is larger than r max , to ensure that all dyons inside the "r max sphere" 
are considered. Since spheres associated with points r near the boundary reach into neighboring 
volumes, one needs to shift the LEGO® sphere periodically. 

4.2 Evaluating the long-range part 

Introducing the structure functions 5(k) = E"=i qje~' krj the long range part can be written 
according to 

v k^o k 

Due to the exponential suppression by e^^^l 1 , it is sufficient to consider momenta 
k < fc max <x\/X. Consequently, one has to sum over &{L 3 jX ) different momenta k and the cost for 
computing the structure functions is ff(V 2 /X 3 ). Since the structure functions 5(k) do not depend 
on r, they have to be computed only once for a given set of dyon positions. 

In a second step the long-range part has to be evaluated at M oc y points r amounting again to 
a computational cost of €?(V 2 /X 3 ). 

4.3 Optimizing the trade-off between the short-range and the long-range part 

The computational cost for the whole algorithm is minimized, when the scaling of the short- 
range and the long-range part are identical, i.e. if 0(VX 3 ) = ff(V 2 /X 3 ). This can easily be achieved 
by choosing A « V 1 ' 6 « \fL, resulting in a total computational cost of ^(V 3 / 2 ). A more detailed 
discussion, in particular, of how to determine the optimal value for X for a specific dyon ensemble, 
we refer to [§]. 

5. Numerical results 

The methods discussed in the previous section were used to evaluate the Polyakov loop corre- 
lator for many different dyon ensembles (cf. [Q] for a detailed discussion). Selected results for 
the quark antiquark free energy as a function of the separation in the case of non-interacting 
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Quark Antiquark Separation dT 



Figure 2: The LEGO® sphere. Figure 3: The quark antiquark free energy as a func- 

tion of the separation d T for T 3 / p = 1 and finite pe- 
riodic volume of length L. 

dyons are shown in Fig. ||. These results correspond to p/T 3 = 1 and four different volumes 
LT e {10,20,30,40}, which is equivalent to n D £ {1.000, 8.000, 27.000, 64.000} (we express all 
dimensionful quantities in units of the temperature T). The curves grow linearly at large quark 
antiquark separations dT and seem to converge to an infinite volume curve with increasing LT . 

To determine this infinite volume curve, we perform a linear % 2 minimizing fit in l/LT for 
each quark antiquark separation and extrapolate to infinite volume, corresponding to l/LT = 
(cf. Fig. |]a). The colored points at l/LT = are analytical results derived in [^|. The numerical 
extrapolations are in agreement to these analytical results within statistical errors. The infinite 
volume free energies (both the numerical Ewald result as well as the analytical result) are shown 
in Fig. ^|b, again demonstrating that by means of Ewald's method one can reliably and efficiently 
determine an infinite volume quark antiquark potential in a long-range dyon ensemble. 



6. Generalization of Ewald's method for arbitrary 1 /r p long-range potentials 

In the previous section we presented numerical results for non-interacting dyons. A more 
realistic model would be to take dyon interactions into account originating from their moduli space 
metric. The problem of non-positive-definiteness of the moduli space metric proposed in can be 
cured by only considering two-dyon, but no three-dyon, four-dyon, etc. interactions fl3|] . One 
then obtains an effective dyon action 



V(k--rjl) 

This effective action also contains "long-range potentials" y, which can be expanded in a power 
series with respect to inverse dyon separations r, 
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Figure 4: a) Quark antiquark free energy as a function of the inverse length of the volume l/LT for different 
separations dT and extrapolations to l/LT = 0. The solid curves indicate the fitting range, whereas the 
dashed lines mark the extrapolations, b) Comparison of the numerical and the analytical infinite volume 
quark antiquark free energy at these. All data was obtained using p /T 3 = 1. 

The individual terms \/r p , p > 1 of this series may be treated with a generalization of Ewald's 
method to arbitrary powers 1 jr p (cf. e.g. Jl4|]), where the short- and long-range parts are given by 

7T 3 / 2 " D / \ / k} \ 

^° ng (r) = — 1 1 qj(p) exp (ik(r - r y )) 4 , (6-4) 



where the charge qj(p) of dyon j may depend on the power p. Again exponential convergence is 
guaranteed, due to 

2 

r(W2) 

fp{x) = Trml s2 ~ Pex p(~ s2 ) ds - (6 - 6) 



2 Z" 30 

8p( - x) = TMi) L ^ lex pH 2 ) d *> ( 6 - 5 ) 

2^;P- 3 



7. Summary and Outlook 

We applied Ewald's method to simulate non-interacting long-range dyon ensembles numeri- 
cally. We have demonstrated that this is an efficient method suited to extract the infinite volume 
quark antiquark free energy. A generalization of Ewald's method offers the possibility to also sim- 
ulate interacting dyons. Such investigations might help to understand the effects and implications 
of the moduli space metric of dyons on the quark antiquark free energy and their relevance regard- 
ing the phenomenon of confinement. Another interesting aspect would be a possible generalization 
to manifestly non-Abelian and typically non-rotationally invariant objects like regular gauge in- 



stantons, merons, or meron pairs, which also seem to yield confinement (cf. e.g. [15, 16, 17]). 
An alternative method to enforce periodic boundary conditions for such non-Abelian long-range 
objects has already been proposed and tested in [p^j. 
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